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Abstract 

We compute the neutral bosonic excitation spectra of the holographic quan- 
tum Hall fluid described by the D3-D7' system. We find that the system is 
stable, gapped, and, in a range of parameters, exhibits a spectrum of low-lying 
excitations very similar to that of a fractional quantum Hall fluid, including a 
magneto-roton excitation. 
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1 Introduction 

Strongly interacting fermions in 2 + 1 dimensions exhibit many interesting phenom- 
ena which are not easily accessible with conventional theoretical tools. Gauge/gravity 
duality offers a new tool which may be better suited for describing and maybe advanc- 
ing our understanding of these phenomena. In [1] , a holographic system exhibiting a 
mass gap at nonzero density and magnetic field was described. This system appears 
to have many of the properties of a fractional quantum Hall fluid. Other holographic 
models of the quantum Hall effect can be found in [2]. 

The fluctuation spectra of quantum Hall fluids have many interesting proper- 
ties, such as gaps for both charged and neutral excitations and fractional charged 
excitations with anyonic statistics. One interesting feature, which has not received 
much attention from holographic studies, is that the lowest energy neutral excita- 
tion is a magneto-roton. Magneto-rotons are collective excitations whose dispersion 
relations have a minimum at some nonzero momentum. Although at this minimum 
the magneto-rotons carry momentum, their group velocity vanishes. Experimentally, 
these are seen in inelastic light scattering experiment [3] and ballistic phonons [1], 
and most recently in [5]. Theoretically, such magneto-roton excitations have been 
analysed using an analogy with the Feynman-Bijl wave-function for the collective 
excitation of the superfluid in [6]jl] the formulation in terms of composite fermions 
in [7], as a Chern-Simons effective theory in [8]. 

In this paper we study the neutral bosonic excitation spectra of the D3-D7' system 
described in [T]. We flnd that the system is stable and gapped. We also flnd that 

^Terminology for the roton is borrowed from studies of the excitations of superfluids, where the 
phonon dispersion has a local minimum at a high momentum. 
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the lowest energy excitation in a certain range of densities displays a magneto-rot on 
minimum at a finite momentum. 

The sectioning of the paper is the following. In Section [2] we present the model 
of [T] in a set of coordinates more suitable for fluctuation analysis. In Section [3] we 
lay out the equations of motion of the fluctuations and the method for solving them. 
In Section H] we present the results. We conclude with a few comments. 

2 Review of the model 

In pp a holographic dual to a (2+l)-dimensional system of strongly-coupled fermions 
in the fundamental representation was proposed. The brane configuration in [1] con- 
sists of a probe D7-brane in the background of — ?■ oo D3-branes. The intersection 
of the D7-brane with the D3-branes is (2+l)-dimensional, such that ^NN = 3 and 
a^ND = 6, and hence the massless spectrum of the D3-D7 open strings consists only 
of the desired fermions (see [5] for earlier consideration of this model, and for its re- 
lationship to three-dimensional QCD, see [lO]). One can then give the fermions mass 
by separating the branes in the transverse direction. 

We decouple these fermions from the closed string modes by focusing on the near- 
horizon limit of a large stack of D3-branes and treating the D7-brane in the probe 
limit. According to the usual gauge/gravity duality, the dynamics of the probe D7- 
brane in the near-horizon D3-brane background captures the physics of these (2+1)- 
dimensional fermions interacting with a strongly-coupled (3+l)-dimensional gauge 
theory. We will review the details of this D3-D7' model and how to obtain stable 
Minkowski (MN) embeddings, which do not enter the horizon, of the D7-brane probe 
at zero temperature. 

2.1 Action 

The D3-brane background has the metric 



The AdS radius is given by = y/AngJV a' = y/Xa', and the temperature is T = 
In addition, the RR four-form Cjf^^^ = "l^- The angles have the following ranges: 
e [0, vr/2], 9i^2 £ [0, tt], and 0i^2 ^ [0, 27r]. In the following, we will only be interested 
in zero temperature, so we set 




(1) 



(2) 



rT = . 



(3) 
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With a view towards the fluctuation analysis, it is useful to introduce the following 
Cartesian coordinates: 

p = rsinip (4) 
R = rcosip . (5) 

We take the D7-brane probe to be wrapped on the two S'^'s and to span the t,x,y, 
and R directions and assume that the embedding fields z and p are only functions 
of R. In order to have a stable configuration, we need to turn on a magnetic flux 
through one of the [Dili]; we choose 

27ra'F02^2 = sin 62 . (6) 

It turns out that with flux only on one of the S'^'s a constant z is a solution to the 
equations of motion; we can thus set it to zero in what follows. The induced metric 
on the D7-brane is 

dsl, = ^ i-dt' + dx' + dy') + ^ (1 + p") dR' + ^dn^;^' + ^dn^^^' . (7) 

Since we are interested in obtaining a quantum Hall description, we will turn on a 
background charge density and a magnetic field: 

27ra'Fr,y = b , 2TTa'Fm = % . (8) 

The D7-brane DBI action reads: 



S = -M I dRR\lP + A^Jl + ^^^/l + p'^- a'o^ , (9) 



where Af = Sn'^T^VsL'^. The CS action is P 

Scs = +ML^ j dR ba'oc{R) , (10) 
where the pull-back of the RR-potential onto the worldvolume reads 

c{R) = arctan (^j^ " ^ sin ^4arctan (^^j j ~ "^oo + ^ sin(4'?/;oo) • (H) 
The asymptotic angle ipoo = hm/^^oo arctan(-^) is related to / by 

f = 4sin2^oo -SsinVoo . (12) 
A rescaling of the coordinates and the temporal component of the gauge field 

R = Ly/bR , p = LVbp , r = Ly/br , = Ly/bdo , (13) 
enables us to write down the full D7-brane action in a ' magnet ic-fieldless' fashion 



S = -Af JdRi^R'^P + A^^l + ^^^l + p - ~a!i + ~a!^c{R)]^ , (14) 
where M = b^^^L^jV and the prime now indicates derivative with respect to R. 
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2.2 Equations of motion 

The a'o equation of motion, integrated once, reads: 



^ 1 



2c{R) 



d 
h 



where 



2R\ P + Af, 



nl2 _ n'-2 



(15) 



(16) 



The constant of integration d was chosen to match the notation of [1] and is related to 
the physical charge density by D = {2Tia' / L){N' /V^)d. After a little bit of an algebra, 
equations f[T^ and f jTBj) can be put in the form 



1 + f-4 

6 



(i 



-\ 2 



+ 4i?4(i + r-4)(/2 + 4|l; 



1 + p 



nl2 



(17) 



1 + p 



The p equation of motion is: 



-16 



^6 



(1 + ^) 



(19) 



We are interested in Minkowski embeddings, for which the D7-brane does not enter 
the horizon. These are characterized by the IR boundary condition p{R = 0) = p^. 
We solve the p equation of motion numerically, shooting from the tip -Rq = 0, using 
a smooth boundary condition p'(0) = 0, toward the UV boundary. As shown in [1], 
the "fermion mass" m and "chiral condensate" c are holographically encoded in the 
UV expansion of 

V' = '0OO + mr^+ - cr^- + . . . , (20) 



where 



A, 



48 



COS^ Ipo 



(21) 



In the rescaled Cartesian coordinates used here, the rescaled mass, m = m/{L\/h), 
can be obtained from the numerically-obtained solution for p{R): 



m 



lim r ^+ sin I arctan 



R—^oo 



(22) 



^The m and c are not really mass and condensate for A-|_ 7^ — 1. 
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As explained in [T], for Minkowski embeddings, the charge density d and the 
magnetic field h are not independent but locked: 

d , ^ TIV , ^ 

- = 2c(0) = — , (23) 

where u/N is the filling fraction per fermion species, which is fixed by the internal 
fiux /. The quantization of the filling fraction is then a direct consequence of the 
quantization of the wrapped fiux. The equation of motion for p f|T9l) fixes ipoo in 
terms of /, which, through equation fl2T]) . also fixes A_|.. Recall from [1] that the 
filling fraction is restricted to lie in the range 

0.6972 < ^ < 0.8045 (24) 

and does not take the usual rational values considered in experimental setups. Par- 
ticular values of the filling fraction for which we will present numerical results are: 

A+ =-1 : ^^ 0.7426 (25) 

A+ = -1 : ^ ^ 0.7082 (26) 

A+ =-| : ^^ 0.6999. (27) 

The dependence of m on po is illustrated in Fig. [T] for three different choices 
of A+@ Note that for fixed (non-rescaled) mass m, lower values of m correspond to 
smaller background magnetic fields b. For a given rn, the stable Minkowski embedding 
is the one with lower po, as was shown by energetic considerations in pQ. We will 
confirm this by showing that the larger po Minkowski embedding has a tachyon in 
its spectrum, whereas for the smaller po it is absent. At 5 = only the unstable 
Minkowski embedding exists. 



3 Fluctuations 

We wish to obtain the spectrum of small fiuctuations around the classical solution 
p{R), ao(-R), and dy{x) = ^x. The energy of the lowest normalizable fluctuation gives 
us the energy gap above the quantum Hall fluid. As before, we write the equations of 
motion and background solutions in the rescaled variables (fT3|) and rescale also the 
fluctuations: 

Sa^ = VbL5d^ (28) 
6p = VbLSp , (29) 

■^We expect there to be also zero-tempcraturc "black hole solutions" reaching the origin of AdS, 
taking any values of the mass [I] . 
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Figure 1: m vs po for MN embeddings. The solid line is the stable MN embedding, 
and the dashed line is the unstable MN embedding. The top, red curve represents 
A+ = —1/2, the middle, blue curve A+ = — 1, and the bottom, black curve A+ = 
-5/4. 



where fi = t,x, y, R. We need to allow for fluctuations of (almost) all worldvolume 
fields: 



z 


= eSz{t,x,y, R) 


(30) 


p 


= p{R) + e6p{t,x,y,R) 


(31) 


at 


= ao{R) + e6dt{t,x,y, R) 


(32) 


ax 


= e5d^{t,x,y,R) 


(33) 




Vb 


(34) 


tty 


= ——x + e6ay{t,x,y,R) 
Ij 




= e6dj^{t,x,y,R) , 


(35) 



with e a small parameter. We do not consider fluctuations in the internal spheres or 
dependence on the internal two two-sphere coordinates. From now on, we will omit 
the bars from the background fields. 

We expand the Lagrangian to second order in e: 

JZ = Co + eCi + e^C2 + ... . (36) 

Here Cq corresponds to the background on-shell Lagrangian and Ci vanishes by the 
background equations of motion. From £2 one can derive the equations of motion for 
fluctuations. 

Not all of the above fluctuations are physical, however. Since we are working in 
the gauge = 0, we must demand that 6dj^ = 0, too. The equation of motion for 
6dj^, after imposing 6dj^ = 0, gives us a constraint which the other fluctuations have 
to obey. 
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We are interested in wavelike solutions for the fluctuation spectrum: 



5z ~ ^-iu)^''h+ik)^'x+ik\;'y ^g-j?^ 

5a^ ~ e-^^i+ik.x+ikyy ^3g^ 

^-iujt+ik^x+ikyy ^ j^gg-j 

where, anticipating the decoupling of 6z fluctuation from the rest, we introduce a 
separate frequency co^^^ and two-momentum k^^^ from the others. All the other fluc- 
tuations couple to one another at nonzero momentum. Rescaling the three-momenta 
as well, 

Co^'^ = ^u:^'^ , ^S = 4^S (40) 



OJ = —uj , = —^k^^y , (41) 



makes the resulting equations of motion dimensionless. 

We also make use of the rotation symmetry in the (x, ?/)-plane and choose the 
vectors k^^^ and k to be aligned along x-axis. From now on, we omit the subscript x 
from the momenta. 

3.1 Equations of motion 

To make sure that all the fluctuations are physical, we write the equations of mo- 
tion in a manifestly gauge invariant way. That is, we introduce a gauge invariant 
combination: 6ex = {k6at + cD^a^.). The equations of motion for 6dt and 6ex are 

|(l + p'^)Me. + |r'-^^ 



df^H = -ikdciydp, [ga'^il + 2f-^)) - ^(1 + P'^)k5e, + ^doP'k^6p 
{^{^^'x ~ kdalf)^ = —iCjdciydf^ (^^o(l + 2f~^)) — ^(1 + p'^)Cj5ex + ^aQp'ujkdp 
The equation for 5ay is 

i5~exdp, (K(l + 2r~')) - + p'')5~ay - {g6~a'y) + ~k' ^ ) ^ J Sdy 

= ikdpdfi [^Pj+ ^/ ao6p - ikN6p , (42) 
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and for 5p we have 

'g 



di 



M - dp,[K~p'] 



1,2 9_ ~/2\ 
4 



I (1 + r-') ~a',p'5~a[ 



+ 47a'Qp'k5e^ + ik I -A^ 



^6 



P «0 + 



(43) 



As mentioned previously, the 5z equation of motion decouples from all the rest: 

[g {f^ + 1) 6z'] = (1 + f-^) (1 + + . (44) 

Finally, the constraint reads 



—g ( 5e^ — k5a[ ) — uH = . 

UJ 



(45) 



In the above equations, we have introduced the following functions which depend only 
on the background: 



1 I ~/2 ~/2 
1 + p -Oq 



(46) 
(47) 



^^14 (1 _^ ^~4) 



2 (/'r^(l + 2f^) - 4fV(l + f^) + 4p^(2 + 3f^)) p (48) 
/^fi°(l + f^) - 2ff'^p\Q + 2/2 + 3f^(4 + f) + 6fS) 



A = 

k = 

N = 
M = 

^^20(1+^-4)2 (^/2+4g 

+ 12/2fV(5 + 8f^ + 3f^) - 8fV(2 + 9f + 2f\2 + if) + f^(2 + 3/')) 
+80f2p^(2 + 3f^ + f^) - 32p^°(6 + 9f^ + 2f^) j . (49) 

In addition, the fluctuation- dependent function H is defined by 
H 



l{i + r%i + p'2)(5a; + 4(1 + f-')~a'oP'5p' 



A 

d3^2 

+16^|^5p + ka'^dp 



A 



(50) 



Looking at the above equations, we notice that after using the constraint, the 
equations of motion are not independent. Indeed, one can see that using the constraint 
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to eliminate 5a[, the equations of motion for 5at and 5ex become identical. So, in fact, 
we are left with only three equations for three unknown functions. Alternatively, one 
can solve four equations for all four functions and impose the constraint only on the 
boundary conditions. 



3.2 Methodology 

We are seeking those solutions which represent normalizable modes. In this parametriza- 
tion, normalizability is translated to the requirement that at the boundary, as -> oo, 

5Sy -> , ^ , R'^~^+Sp . (51) 

At the tip, the physical boundary conditions are that all derivatives with respect to 
R vanish: 

Sp'iO) = 6e'M = K(0) = . (52) 

The vanishing of Sp'{0) means that the embedding is smooth, and the boundary 
conditions for the gauge fluctuations ensure that there are no sources at the tip. 

The equations of motion (other than for 6z) cannot be completely decoupled|f| 
which makes our goal of finding solutions which obey normalizable boundary condi- 
tions at infinity much more difficult. Operationally, we solve the equations by choos- 
ing boundary conditions at the tip and shooting out to infinity. To find the desired 
normalizable solutions, we must systematically vary the choice of tip boundary con- 
ditions until we find those for which all the coupled fiuctuations will simultaneously 
vanish at infinity. 

To accomplish this task, one can utilize a known determinant method [T^[T5] . We 
have specified three boundary conditions at the tip for the derivatives of the fiuctua- 
tions. Since, we have three independent functions to solve for, and they satisfy ordi- 
nary, linear, second-order homogeneous differential equations, we should choose three 
linearly-independent boundary conditions for the fluctuations themselves. Following 
the method of [12], we choose a set of linearly-independent boundary conditions, e.g., 

(<5p(0), 5e,(0), 5~ay{0)) = {(1, 0, 0), (0, 1, 1), (0, 1, -1)} , (53) 

which span the space of all possible boundary conditions at the tip. We then solve 
the equations of motion numerically for a fixed (a), k) for each group and compute 
the following determinant: 



R' 



det 



^+-'6p' 


R-^+-^6p" 










H 







(54) 



*At k ^ 0, Sp decouples from the rest. 
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where the Roman index corresponds to the given group of initial conditions in (153|1 . 
For a given k, we then scan through u (a) can also take complex values) and calculate 
the determinant at each step. When the determinant is zero, there is a linear com- 
bination of the initial conditions giving the desired normalizable solution for which 
all fluctuations vanish at infinity. Said differently, if for the set (a), k) there is such a 
solution, it means that at (a), k) there is a pole in the Green's functions [T2] . 



4 Numerical results 

We first present the numerical analysis for the case k = 0. For this case, in addition to 
the generic decoupling of 6z, the {6ex,Sdy) system and the {6p,6dt) system decouple 
from each other as well. The energy of the normalizable excitation represents the 
gap to the next energy state. At zero k, we can identify {6er^,6dy) as the vector 
excitation and (5p, 6dt) as the scalar excitation, but since they mix at nonzero k, 
such a distinction is not generically possible. Nevertheless, we can continue to refer 
to these modes by their k = labels. 
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Figure 2: The mass (energy gap) at A; = in units of the magnetic length ^ for 
the first few excitations. In the left panel we show the excitations in the vector 
sector, middle panel shows the excitations in the scalar sector, and the right panel the 
excitations in the scalar Sz. All figures are generated at A+ = —1 and m = —2.185. 



Fig. [2] shows the first few excitations of the various fields. First, we do not find 
any tachyons, implying that the embedding is stable. We also see that the embedding 
describes a gapped state in the neutral sector which means that it is an incompressible 
fluid. The spectrum of the excitations is approximately evenly spaced with a scale set 
by (g2 ^)i/4 ; where B is the physical magnetic field. Note that in the vector sector, 
the excitations come in doublets except for the lowest one. 

For most excitations, u grows when the magnetic field increases (at fixed filling 
fraction). However, the energy of the lowest scalar excitation decreases at large 
enough magnetic field. To illustrate this, we plot in Fig. |3] the energy, in units of 
the magnetic length of the lowest scalar mode and the lowest vector mode. At 
small magnetic field, the vector has lower mass, but as the magnetic field grows, the 
scalar mass drops down and eventually vanishes. This occurs for large magnetic fields 
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Figure 3: The mass (energy gap) squared at = in units of the magnetic length, of 
the two lowest energy excitations at zero momentum are depicted in the left panel. 
The red (blue) line represents the scalar (vector) fluctuation. Dashed lines are the 
same fluctuations in the unstable MN embedding. Notice that the scalar fluctuation 
is tachyonic in the unstable MN branch. In the right panel, the energy gap at A; = 
as a function of the magnetic field (or density). The vector fluctuation increases in 
the whole range of the magnetic field. Both figures are generated at A+ = — 1. 

close to where the branch of the stable MN solution meets the branch of the unstable 
solution. Going onto the unstable branch, we see the tachyon explicitly as the scalar 
embedding fluctuation. 



(h 




Figure 4: In the left panel the dispersion relation of the vector fluctuations (blue) and 
of the scalar fluctuations (red, the first and the fourth curve from the top) at m = —3 
and A_|_ = — 1. In the right panel the dispersion of 5z at m = —3 and A_|_ = — 1 is 
shown. 

We now turn to the analysis of the dispersion relations of these excitations, starting 
with the higher mass fluctuations. Their dispersions are plotted in Fig. H] and are of 
the form 

u = \lcjl + clk^ , (55) 

where is the speed of sound for these fluctuations. Interestingly, we find that for 
sufficiently large momentum some of the lines cross each other, which means that 
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Figure 5: The dispersion relation of the two lowest energy excitations for fh = 
—2.15,-2.185, and —2.25 from left to right, at A_|_ = —1. The lowest excitation 
displays a roton minimum u = at nonzero k = k^. 



there is some degeneracy. We also find that as the magnetic field is increased (for 
fixed filling fraction), Cg increases towards one. 

We now focus on the more interesting case of the two lowest energy modes, one 
of scalar origin and one of vector origin. As the magnetic field is increased (for 
fixed filling fraction), the lowest energy scalar mode decreases in energy, while the 
vector mode increases slightly. As depicted in Fig. El at some critical magnetic field 
the lowest energy mode becomes a roton; this is, it develops a minimum at nonzero 
momentum. The dispersion assumes the characteristic rotonic form: 

u = ^JcJl + cl{k-~Kf , (56) 

where the star corresponds to the values at the lowest point of the dispersion curve 
for given fh. Although at this minimum the roton carries a nonzero momentum /c*, 
its group velocity is zero. 

As the magnetic field is increased further, the roton minimum a)* initially gets 
deeper, then k^ begins decreasing, until finally the minimum disappears as K — ?• 0. 
We depict this behavior in Fig. [6] in units of the magnetic length. In this range of fh 
where the roton exists, it is always the lowest energy excitation of the system. Notice 
that the roton only exists in the range of parameters where there is also a neighboring 
excitation. This suggests that the near degeneracy is correlated with the existence of 
the roton. 

This behavior is seen for all filling fractions, and the region of parameter space 
where the roton exists is around the region of parameter space where the energy 
of the two lowest excitations become close. The existence of two very close energy 
excitations has been recently observed in experiments [2] and has received some 
theoretical explanation from hydrodynamics [15] and from the composite fermion 
picture [T6|fT7]. 

In general, we find that as the filling fraction decreases, the minimum energy 
of the excitation decreases as well. In addition, the value of fh at which the roton 
appears also decreases with decreasing filling fraction. To illustrate this, we plot in 
Fig. [7] the dispersion curves for various values of v/N , with fh chosen for each one 
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Figure 6: The location of the roton minimum as a function of 1/fh? is shown in the 
left panel. The right panel displays the minimum of the roton energy as a function 
of 1/fh?. Both figures are for = —1. 
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Figure 7: Dispersion relations for the two lowest energy states when their energy 
is degenerate at zero momentum: (left) = —1/2,-^ ^ 0.7426, (middle) A_|_ = 
-1, ^ ^ 0.7082, and (right) A+ = -5/4, |r ^ 0.6999. Notice that the minimum 
energy cD* decreases with decreasing filling fraction. 

such that the two lowest excitations are degenerate at zero k. We see that goes 
down with decreasing v/N but does not reach zero in the range of (12^ of possible 
filling fractions. 

5 Discussion 

In this paper we continued to investigate the holographic description of a quantum 
Hall fluid first described in [1]. This D3-D7' model was shown previously to have 
the zero longitudinal conductivity and the correct Hall conductivity appropriate for 
a quantum Hall fluid. Here we analyzed the neutral sector and showed that it is 
gapped as well, leading to an incompressible fluid. We also showed that the lowest 
neutral excitation has a dip in its dispersion relation leading to its identification as a 
magneto-roton. 

Another intriguing similarity with observed quantum Hall systems is the existence 
of an excitation with ordinary dispersion slightly above the rotonic excitation (as 
shown in Fig. E]). At small momentum, the non-rotonic excitation was observed 
in [H] to be very close to the rotonic one. In our model the two are close only in a 
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range of densities; it is unclear in real materials what role the density plays in this 
picture. We also note that it is believed that the non-rotonic excitation is a two-roton 
bound state [l7j, while in the D3-D7' model, it is just one of the mesons. 

While we found the model agrees with many expected properties of the quantum 
Hall fluid, there are also some significant differences. In terms of the fluctuation 
spectrum, although we have reproduced a magneto-roton, it does not completely 
match with those observed in experimental quantum Hall systems. 

In this model there is a large hierarchy between the lowest energy neutral ex- 
citation, which is of the order of A/5((7yjjA^)~^/^, and the energy scale of charged 
excitations, which is expected to be of the order of ^/B{gYM^y^^- This is rather dif- 
ferent from real quantum Hall systems where there is no such parametric gap. Thus 
in the holographic model, where Qym^ is large, the lowest neutral excitations are 
very deeply bound quasiparticle-quasihole bound states. This difference in energy is 
geometrically clear from the point of view of the bulk; quasiparticles are expected 
to be formed by long strings stretched from the tip of the D7-brane to the horizon 
(probably bound to some wrapped five branes), and the lowest neutral excitations 
are short strings representing fluctuations of the D7-brane. 

The magneto-roton we found seems to exist only for a range of magnetic fields (or 
densities), and we found just one minimum in the dispersion relation for any filling 
fraction. This is in contrast to real materials where the number of roton minima 
depends on the filling fraction and, as far as we know, is not sensitive to the density. 

Although our results support the proposal of [1] that the D3-D7' system holo- 
graphically describes a quantum Hall fluid, it is an imperfect model. This is not 
entirely surprising, considering that holographic models necessarily have various un- 
physical properties, for example large A^, which are undesirable from the point of view 
of applications to physical systems. However, because the quantum Hall effect is a 
generic feature of (2+l)-dimensional fermion systems with broken parity and exhibits 
many universal properties, we are ultimately optimistic about our chances to describe 
it holographically. 
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